Linear Magnetic Response of Disordered Metallic Rings: Large Contribution from 

Forward Scattering Interactions 
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We calculate the effect of electron-electron interactions involving vanishing momentum transfer 
(forward scattering) on the orbital linear magnetic response of disordered metal rings pierced by a 
magnetic flux (j). Using the bulk value of the Landau parameter _fo for copper, we find that in the 
experiment by Levy et al. [Phys. Rev. Lett. 64, 2074 (1990)] the forward scattering contribution 
to the linear magnetic response is larger than the corresponding contribution from large momentum 
transfers considered by Ambegaokar and Eckern [Phys. Rev. Lett. 65, 381 (1990)]. However, 
outside the regime of validity of linear response and to first order in the effective screened interaction 
the persistent current is dominated by scattering processes involving large momentum transfers. 
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PACS numbers: 73.23.Ra, 74.20.Fg 

I. INTRODUCTION 

More than a decade ago the measurement by Levy et 
al} of persistent currents in mesoscopic normal metal 
rings pierced by an Aharonov-Bohm flux (p has triggered 
a lot of theoretical activitj*2i^. Yet, up until now a truely 
convincing and generally accepted theoretical explana- 
tion of the surprisingly large persistent currents observed 
in Ref.^ and in subsequent experiments^'^ has not been 
found. It has become clear, however, that this effect 
cannot be explained within a model of non-interacting 
electrons. Ambegaokar and Eckern (AE)^ were the first 
to examine the effect of electron-electron interactions on 
mesoscopic persistent currents: they realized that, to first 
order in the screened Coulomb-interaction, the dominant 
contribution to the disorder averaged persistent current 
can be obtained from the two diagrams shown in Fig. 
n representing a special correction riAE(0) to the dis- 
order averaged thermodynamic potential which depends 
strongly on the Aharonov-Bohm flux (j). Here the over- 
line denotes averaging over the disorder. Given the grand 
canonical potential 11(0), the corresponding persistent 
current I{(j)) can be obtained from the thermodynamic 
relation 



(1) 



In a bulk metal at high densities the bare Coulomb- 
interaction Vb(q) = 47re^/q^ is strongly screened. A sim- 
ple way to take the screening into account diagrammati- 
cally is the random-phase approximation (RPA) . Follow- 
ing this procedure, AE approximated the effective inter- 
action (in the imaginary frequency formalism) as follows 
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FIG. 1: Feynman diagrams representing the flux-dependent 
part of the grand canonical potential to first order in the 
screened interaction, (a) Hartree diagram; (b) Fock dia- 
gram. Solid arrows represent non-interacting disorder av- 
eraged Green functions and thick wavy lines represent the 
effective density-density interaction. The Cooperon (shaded 
symbol) is defined in Fig. |5| 
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FIG. 2: Diagrammatic definitions of the Cooperon (C) and 
the Diffuson (D). A dashed line represents the covariance of 
the impurity potential. 
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For momentum transfers |q| small compared with the in- 
verse elastic mean free path £~^, and for frequency trans- 
fers \uj\ small compared with the inverse elastic lifetime 



T ^ the disorder averaged polarization is given by 
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where Dq is the diffusion coefficient and vq is the average 
density of states at the Fermi energy (per spin) in the 
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absence of interactions. Note that = (AoV)"""^, where 
V is the volume of the system and Aq is the average 
level spacing (per spin) at the Fermi energy. It turns out 
that both diagrams in Fig. ^^re dominated by momen- 
tum transfers of the order of the Fermi momentum kp, 
which for a metallic system is large compared with £^^. 
Eqs. |2Jl and Q are therefore not suitable for a quanti- 
tatively accurate calculation of persistent currents. To 
make some progress analytically, AE estimated the con- 
tribution from the diagrams in Fig. ^ by replacing the 
effective interaction by a constant 

FRPA(k - k',tiu) ^ (FRPA(kf^ - k'p,tO)) = V, (4) 

where (...) denotes the Fermi surface average over kj? 
and k^. For simplicity, it is assumed that the ring is 
quasi one-dimensional, with transverse thickness L± in 
the range kp"^ <C L±_ i L, where L is the circumfer- 
ence of the ring. Then diffusive motion is only possible 
along the circumference. At temperature T = the re- 
sulting average persistent current can be written as& 



^4^^^sin(4^fc0/(/.o) 
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where (f>Q = hc/e is the flux quantum^ and the Fourier 
coefficients of the current are 

c 16Ac 
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Here Ec ~ HDo/L^ is the Thouless energy and 7 = 
T/Ec ^ 1, where at zero temperature F = Aq/tt is 
the cutoff energy that regularizes the singularity in the 
Cooperon in a finite system^, see Eqs. (|14|l and H2U|) 
below. The coupling constant Ac = i'qV can be iden- 
tified with the dimensionless effective interaction in the 
Cooper channel to first order in perturbation theory. AE 
estimated Ac ~ 0.3, assuming that the validity of the 
RPA can be extended to momentum transfers of the or- 
der of fci?- However, higher order ladder diagrams in the 
Cooper channel strongly reduce the effective interaction, 
so that Ac ~ 0.06 is a more realistic estimated for the 
Cu-rings in the experiment. 

In real space Eq. Q amounts to replacing the electron- 
electron interaction by a local effective density-density 
interaction. 



T/eff(r-r') ^ VS{r~r') 
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More precisely, this replacement means that for distances 
|r — r'l larger than £, the interaction is effectively local. 
In a recent letter Schechter, Oreg, Imry, and LevinsorkiS 
pointed out that a different type of effective interaction 
can possibly lead to a much larger persistent current. 
Specifically, they used the BCS model to calculate the 
leading interaction correction to the orbital linear mag- 
netic response and founc&ifi 



= --2 327rABcs£^cln 



6=0 



where Abcs < is the attractive dimensionless interac- 
tion in the BCS model, and the coherence energy E^o 
is the smaller energy oi h/r and the Debye energy fuoo. 
Eq. ||SJ| should be compared with the corresponding re- 
sult for the local interaction model used by AE, which 
implies according to Eqs. (0) and jnj. 
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where we have used F — Aq/tt and retained only the 
leading logarithmic order. Note that the logarithm is 
due to the slow decay (ex fc^^) of the Fourier coeffi- 

— AE 

cients inkl^^ /(j)o of dl /d(p, so that all coefficients 
with k < 1/^7 contribute to the linear response. For 
Eco ^ Ec the linear magnetic response in the BCS model 
is parametrically larger than the linear response in the lo- 
cal interaction model. Whether or not this remains true 
beyond the linear response has not been clarified. Note 
also that in the BCS model the linear magnetic response 
is diamagnetic because the effective interaction is attrac- 
tive (Abcs < 0), whereas the linear response in the local 
interaction model is paramagnetic, corresponding to a 
repulsive effective interaction (Ac > 0). 



II. MAGNETIC RESPONSE DUE TO 
FORWARD SCATTERING 

An interesting observation made by the authors of 
Refiiii is that an effective interaction different from the 
local interaction used by AE can lead to a much larger 
persistent current, at least for sufficiently small flux 0, 
where it is allowed to calculate the current from the lin- 
ear response. Given the rather crude approximations in 
the microscopic derivation of the local interaction model, 
it seems worth while to explore the magnetic response for 
other types of effective interactions. A possibility which 
so far has not been thoroughly analyzed is an interaction 
which is dominated by small momentum transfers. Note 
that the assumption that only forward scattering pro- 
cesses (corresponding to vanishing momentum transfer) 
have to be taken into account for a consistent descrip- 
tion of the low-energy and long-wavelength properties of 
normal metals lies at the heart of the Landau's Fermi 
liquid theory. The Landau model is in a sense the oppo- 
site extreme of the local interaction model, because the 
effective interaction in the Landau model is proportional 
to a Kronecker-delta in momentum space. 



V^off(q, «t^) 5q,o/o 
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where the Landau parameter /o can be determined from 
experiments. In fact, the dimensionless Landau param- 
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Eq = 2voJo can be written as Fq = "i:^ - 1 



Bo mo 



where B is the bulk modulus, m» is the effective mass, 
and Bq and mg are the corresponding quantities in the 
absence of interactions. Inserting the known bulk values 
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where 



FIG. 3: Feynman diagrams that dominate the flux-dependent 
part of the grand canonical potential if the effective inter- 
action involves only momentum transfers smaller than the 
inverse elastic mean free path, |q| < . For vanishing mo- 
mentum transfer the Hartree diagram (a) dominates the linear 
magnetic response, whereas outside the regime of validity of 
linear response the sum of the three Fock diagrams (b)-(d) 
has the same order of magnitude as the Hartree diagram (a). 
Note that the Diffuson (shaded box, see Fig.|5J renormalizes 
only the density vertex in the Fock diagrams; the interaction 
in the Hartree diagram does not transfer any energy and hence 
cannot be renormalized by singular Diffuson corrections. 
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Essentially we have used the equation of continuity to 
replace the charge vertices in Fig. |31by current vertices, 
which cannot be renormalized by singular diffusion cor- 
rections. The fact that a gauge transformation replacing 
charge vertices by current vertices can be used to avoid 
the explicit calculation of vertex corrections has also been 
employed in Ref.^^ to calculate the zero bias anomaly in 
the tunneling density of states of two-dimensional disor- 
dered electrons interacting with Coulomb forces. 

The evaluation of the contribution of the three Fock di- 
agrams in Fig.|31to the persistent current is now straight- 
forward. Note that for a thin ring with L± <C ^ <C i the 
q-summation is one-dimensional, with quantized wave- 
vectors 27r(n -I- 2(t)/<j>Q)/L, n = 0, ±1, ±2, . . .. Then we 
obtain for the fc-th Fourier component of the average cur- 
rent due to the Fock diagrams (b)-(d) in Fig. |31 for the 
Landau modeli^, 



for Cu^^ TO,/m « 1.3 and B/Bq « 2.1, we find Fq « 1.7, 
which is a factor of 30 larger than the corresponding es- 
timate Ac « 0.06 in the local interaction model. Note 
that in real space Eq. (|10|) corresponds to a constant ef- 
fective interaction, proportional to the inverse volume of 
the system 
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Given an effective interaction of the form H10|l , the dom- 
inant flux-dependent contributions to the average poten- 
tial to first order in the interaction are shown in Fig. 
13 The Fock diagrams (b)-(d) have been discussed pre- 
viously in Refsii2*i^; as first pointed out by Beal-Monod 
and Montambaux^"^, to leading order in the small param- 
eter {kp£)~^, the three Fock diagrams in Fig. |2| (b)-(d) 
cancel, so that a direct evaluation of the sum of these di- 
agrams is rather difficult. To calculate the leading contri- 
bution of these diagrams, we note that the fermion loops 
in Fig. O (b)-(d) can be identified with contributions to 
the disorder averaged polarization, which for general fre- 
quencies and small wavevectors can be written asi^ 



no(q,it^) = 21^0 



D{iuj)q^ 
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where D{iuj) is a generalized frequency-dependent Diffu- 
sion coefficient. The crucial observation is now that the 
sum of the three Fock diagrams in Fig. |3| (b)-(d) corre- 
sponds to the usual weak localization correction to the 
average conductancei^, 



^L,Fock J.-1/0 

'-k °^ v" 



(15) 



Due to the extra factor of inverse volume, this contri- 
bution is, for experimentally relevant parameters^, neg- 
ligible compared with corresponding result in the local 
interaction model given in Eq. 

The Hartree diagram in the Landau model is more in- 
teresting. The fact that the diagram with two Cooperons 
shown in Fig.|3|(a) dominates the persistent current due 
to electron-electron interactions with momentum trans- 
fers |q| < has already been pointed out in Ref.^^. A 
similar diagram with two Cooperons (but without inter- 
action line) dominates the fluctuations of the number of 
energy levels in a fixed energy window centered at the 
Fermi energyiS. Using the approximate relation 



J ^('^^)' 



(16) 
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between the persistent current In{4>) at constant par- 
ticle number and the fluctuation [5N)'^ of the particle 
number at constant chemical potential /i, several au- 
thors have realizedii2i22i2i that without interactions the 
two-Cooperon diagram determines the average persistent 
current in a canonical ensemble. Note that the Hartree 
diagram in Fig. |31 (a) does not contain any vertex cor- 
rections analogous to the diffusion corrections of the ver- 
tices in the Fock diagams (b)-(d). This is due to the 
fact that the interaction line in the Hartree process does 
not transfer any energy. Hence, the two Green functions 
attached to the vertex of a Hartree interaction are either 
both retarded or both advanced, so that it is impossible 
to attach a singular Diffuson to the vertex. 

For the Landau model the Hartree diagram in Fig. |2| 
(a) yields at finite temperature T the following correction 
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to the disorder averaged grand canonical potential, 
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x^[Go(k',iL;)„0]^Go(-k' + q,it:^„) . 
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Here r = i/vp is the elastic lifetime, (D„ = 27r(n 
are fermionic Matsubara frequencies, and 



Go(k, iCbn) 



1 



^ + M + «|:signi:^„ 



(17) 
i)T 

(18) 



is the disorder averaged non-interacting Matsubara 
Green function. Since the Cooperons (i.e. the second 
line) in Eq. p7|) are only singular for |q| < , and be- 
cause the k- and k'- sums are dominated by momenta of 
the order of the Fermi momentum, we may approximate 
Go(-k-t-q, ia)„') « Go(-k, ia)„') and Go(-k' + q, iCbn) ~ 
Go(— k',iJ;„) in Eq. ((T7|l . The product of the last 
two lines of Eq. (|17|) gives then rise to a factor of 
[(27r/Ao)(r/?i)^]^, so that we obtain 



-:^L,Hartrcc - > 
il {(j)} 
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with the dimensionless coefficient 



0<ui„-,<h/T q 



(19) 
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Note that the small energy scale Aq of Eq. 1)21(1 has disap- 
peared on the right-hand side of Eq. , and is replaced 
by the much larger Thouless energy Ec- Due to the faster 
decay of the Fourier components (|15|) of the Fock contri- 
bution in the Landau model, the linear response due to 
the Fock diagrams shown in Fig. 13 (b)-(d) is a factor of 
y^T/Ec smaller than the corresponding Hartree contribu- 
tion. Interestingly, the anomalously large linear magnetic 
response in the BCS model given in Eq. © is also dom- 
inated by the Hartree processiS. Thus, the importance 
of Hartree interactions for persistent currents is some ex- 
tent independent of a specific model for the interaction. 
For the Cu-rings used in the experimenlji we estimate 
Fo ~ 1.7, Ac ~ 0.06 and Ec/Aq w 25; with these values 
the linear magnetic response due to forward scattering is 
more than four times larger than the linear response in 
the local interaction model considered by AE&. To take 
both contributions into account one should parameterize 
the total effective interaction as 

Feff(r-r')=FJ(r-r') + ^^^, (24) 



V 



which in momentum space amounts to 



V^eff(q) 



V forq = 
fo forq^O 



(25) 



Using the estimate for bulk Cu given above, Fq = 
2vofo w 1.7 and Ac = v^V w 0.06, we find /o/F w 14, 
which supports our assumption that there is indeed a 
strong enhancement of the effective interaction in the for- 
ward scattering channel. 



where lo^ — 2TrmT are bosonic Matsubara frequencies. 
Assuming again a thin ring with L± <^ £ <^ L, we find 
in the limit T ^ for the Fourier components of the 
persistent current. 
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Comparing this expression with the corresponding result 
© of the local interaction model, we see that in the 
Landau model the Fourier components /; 



L,H3-rtree 



are in- 



dependent of k as long as k < 1/y^. Therefore the linear 
magnetic response is de termined by all Fourier compo- 
nents up to fc < Ec/T, 
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where we have assumed that 7 = F/E'c <C 1, so that 

eVi 1 



fc=i 
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III. CONCLUSIONS 

In summary, we have shown that a forward scatter- 
ing excess interaction, involving a vanishing momentum 
transfer, yields the dominant contribution to the linear 
magnetic response of mesoscopic metal rings for experi- 
mentally relevant parameters!^. On the other hand, out- 
side the linear response regime the persistent current is 
dominated by the term V involving large momentum 
transfers, at least if we use the bulk estimates for V 
and Fq for Cu. However, for a mesoscopic disordered 
Cu-ring it is not obvious that the bulk estimates are reli- 
able. Note also that in the bulk the normal Fermi liquid 
is stable as long as Fq > —1, so that for > > — 1 
the linear magnetic response in the normal state can be 
diamagnetic, in spite of the fact that the effective cou- 
pling Ac in the Cooper channel is positive. Moreover, in 
the vicinity of an s-wave Pomeranchuk instability^^'^'^, 
where Fq < and |1 -I- -FqI ^ 1, one should replace Fq 
by ^o/(l + ^o)- In this case we predict a strongly en- 
hanced diamagnetic linear response. In fact, the forward 
scattering channel might then dominate the persistent 
current even beyond the linear order in the flux 0. To 
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clarify this point, a better microscopic theory of the ef- 
fective electron-electron interaction in mesoscopic disor- 
dered metals is necessary. In particular, a microscopic 
theory should properly treat the problem of screening in 
a finite system and incorporate the breakdown of Fermi 



liquid theory in quasi one-dimensional disordered metals 
at sufficiently low temperatures^!. 

We thank M. Schcchtcr for his clarifying remarks con- 
cerning Ref^ and for his comments on this manuscript. 
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